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ABSTRACT 

The efficiency of star formation, defined as the ratio of the stellar to total (gas and stellar) mass, is 
observed to vary from a few percent in regions of dispersed star formation to about a third in cluster- 
forming cores. This difference may reflect the relative importance of magnetic fields and turbulence in 
controlling star formation. We investigate the interplay between supersonic turbulence and magnetic 
fields using numerical simulations, in a sheet-like geometry. The geometry allows for an accurate and 
expedient treatment of ambipolar diffusion, a key ingredient for star formation. We demonstrate that 
star formation with an efficiency of a few percent can occur over several gravitational collapse times 
in moderately magnetically subcritical clouds that are supersonically turbulent. In turbulent clouds 
that are marginally magnetically supercritical, the star formation efficiency is higher, but can still be 
consistent with the values inferred for nearby embedded clusters. A phenomenological prescription 
for protostellar outflow is included in our model to stop mass accretion after a star has obtained a 
given mass and to disperse away the remaining core material. Within a reasonable range of strength, 
the outflow does not affect the efficiency of star formation much and contributes little to turbulence 
replenishment in subcritical and marginally supercritical clouds. If not regulated by magnetic fields 
at all, star formation in a multi- Jeans mass cloud endowed with a strong initial turbulence proceeds 
rapidly, with the majority of cloud mass converted into stars in a gravitational collapse time. The 
efficiency is formally higher than the values inferred for nearby cluster-forming cores, indicating that 
magnetic fields are dynamically important even for cluster formation. 

In turbulent, magnetically subcritical clouds, the turbulence accelerates star formation by reducing 
the time for dense core formation. The dense cores produced are predominantly quiescent, with sub- 
sonic internal motions. These cores tend to be moderately supercritical, and thus remain magnetically 
supported to a large extent. They contain a small fraction of the cloud mass, and have lifetimes long 
compared with their local gravitational collapse time. Some of the cores collapse to form stars, while 
others disperse away without star formation. All these factors, as well as core-outflow interaction, 
contribute to the low efficiency of the star formation in these clouds of dispersed star formation. 
Subject headings: ISM: clouds — ISM: magnetic fields — MHD — stars: formation — turbulence 



1. INTRODUCTION 

1.1. Motivation and Previous Work 

Stars are formed out of molecular clouds that are both 
turbulent and magnetized in the present-day Galaxy. 
The relative importance of magnetic fields and turbu- 
lence in star formation is a matter of debate. One school 
of thought is that star formation is primarily controlled 
by turbulence (Larson 1981), with magnetic fields play- 
ing a minor, if any, role (Mac Low & Klessen 2004). In 
this picture, supersonic turbulence produces dense cores 
in shocks (e.g., Padoan et al. 2001; Gammie et al. 2003; 
Li et al. 2004), some of which can become self-gravitating 
and collapse to form stars. The shock- produced dense 
cores tend to evolve dynamically (Ballesteros-Paredes, 
Klessen & Vazquez-Semadeni 2003), unless the postshock 
region happens to contain a mass close to the local ther- 
mal Jeans mass. The dynamic motions are in contrast 
with the observation that the starless cores of nearby 
dark clouds tend to be quiescent, with subsonic nonther- 
mal velocity dispersions (e.g., Jijina, Myers & Adams 
1999). In the best candidates of prestellar collapse (such 



as L1544), the infall speeds derived from line- profile mod- 
eling are typically subsonic as well (Tafalla et al. 1998; 
Lee, Myers & Tafalla 1999). Furthermore, a strongly 
turbulent, weakly magnetized cloud that contains many 
thermal Jeans masses, if not driven constantly on a scale 
smaller than the Jeans length (Klessen, Heitsh & Mac 
Low 2000), would collapse in a turbulence crossing time, 
forming stars with a high efficiency (e.g.. Bate, Bonnell 
& Bromm 2003). The high efficiency may or may not 
be a problem for localized regions of cluster formation 
(Lada & Lada 2003); it is well above the few percent 
level that is typically observed in molecular clouds as a 
whole (Evans 1999). The problem is particularly severe 
for regions of dispersed low-mass star formation, such as 
the Taurus molecular clouds. 

Inefficient star formation is characteristic of the stan- 
dard scenario of isolated low-mass star formation (Shu, 
Adams & Lizano 1987). In this picture, a strong, ordered 
magnetic field is postulated to prevent the bulk of the 
cloud material from prompt collapse. Stars are formed in 
dense pockets of the magnetically supported, subcritical 
background, where the magnetic support has weakened 
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through ambipolar diffusion (Nakano f984). Many au- 
thors have done detailed ealculations of cloud evolution 
based on this scenario (see Mouschovias & Ciolek 1999 
and references therein), assuming for numerical conve- 
nience that the clouds are non-turbulent. Zweibel (2002) 
and Fatuzzo & Adams (2002) have shown analytically 
that turbulence can speed up ambipolar diffusion in the 
relatively low density regions where turbulent motions 
are dynamically important. 

Whether the bulk cloud material is indeed magneti- 
cally subcritical remains to be determined observation- 
ally. Available Zeeman measurements (Crutcher 1999) 
indicate that the magnetic flux-to-mass ratios in molec- 
ular clouds are not far from the critical value (27r(j^/^; 
Nakano & Nakamura 1978) for magnetic support, al- 
though the exact values are uncertain, because of un- 
certain projection effects (Shu et al. 1999). Shu et 
al. advanced an argument for the molecular clouds in 
the present-day Galaxy being close to magnetically crit- 
ical: those far from the criticality would have either ex- 
panded to become diffuse HI clouds or collapsed to be- 
come the interiors of stars long time ago. In perhaps the 
best studied starless core, L1544, the line-of-sight field 
strength was measured at ^ 11 /xG (Crutcher & Troland 
2000), corresponding to a deprojected flux-to-mass ra- 
tio of about half the critical value, in agreement with 
the prediction of models of core formation driven by am- 
bipolar diffusion (Basu & Mouschovias 1994; Nakamura 
& Li 2003). The agreement lends support to the picture 
of low-mass star formation involving ambipolar diffusion. 
A potential problem of this picture is that it takes an un- 
acceptably long time for ambipolar diffusion to produce 
dense cores directly out of the low-density background 
regions that are well ionized by interstellar UV photons 
(McKee 1989; Myers & Khersonsky 1995). This problem 
can be resolved by strong turbulence, which accelerates 
the ambipolar diffusion-regulated star formation through 
shock compression. 

Detailed investigations of the combined effects of a 
strong turbulence and magnetic field on star formation 
in the presence of ambipolar diffusion require numer- 
ical simulation. In a recent work, we have demon- 
strated, using MHD simulations in a 2D sheet-like ge- 
ometry, that supersonic turbulence and ambipolar diffu- 
sion can work together to convert a small fraction of the 
mass of a moderately magnetically subcritical cloud into 
dense collapsed objects over several turbulence crossing 
times (Li & Nakamura 2004; Paper I hereafter). The 
magnetic fields prevent the turbulence from compress- 
ing most the cloud mass into collapse in one crossing 
time, which would happen if the cloud is weakly mag- 
netized. The turbulence, on the other hand, acceler- 
ates the reduction of magnetic flux in shocked regions. 
In this paper, we will continue our effort in develop- 
ing this "turbulence-accelerated, magnetically regulated" 
scenario of low-mass star formation, which appears par- 
ticularly promising for the relatively inefficient, dispersed 
mode of star formation observed in nearby regions like 
the Taurus molecular clouds. This mode is difhcult to 
accommodate in the purely turbulent picture of star for- 
mation. 

1.2. Protostellar Outflows and Stellar Mass 
Determination 



In Paper I, we took a first step towards determining 
the efficiency of star formation, defined as the ratio of the 
stellar to total mass of stars and gas (Lada & Lada 2003). 
We computed the fraction of cloud mass that has become 
dense (10 times above the average column density) and 
magnetically supercritical. This dense supercritical ma- 
terial can be loosely identified with the dense cores of 
molecular clouds, which are known to be intimately tied 
to low-mass star formation. There is, however, evidence 
that only a fraction of the core material eventually ends 
up inside stars. Onishietal. (2002), for example, derived 
an average virial mass of about 5Mq for the H^'^CO"'" 
cores in Taurus, which is an order of magnitude higher 
than the mass of a typical star formed in that region 
(Kenyon & Hartmann 1995; Palla k Stahler 2002). It is 
unclear how the stellar mass is limited to only a small 
fraction of the core mass. A widely discussed possibil- 
ity is that the mass accretion is stopped by a powerful 
protostellar outflow (Shu et al. 1987; Matzner & McKee 
2000), particularly if the core is surrounded by a mag- 
netically subcritical envelope. The subcritical envelope 
would reduce the mass accretion rate at late times, mak- 
ing it easier for the outfiow to reverse the infall (Shu et 
al. 2004). The infall-outflow interaction appears to be 
observed in B5 (Velusamy & Langer 1998), among oth- 
ers. The details of the interaction remain to be worked 
out. 

In the absence of a detailed theory of how protostellar 
mass accretion is stopped, we will turn to observation 
for guidance. In low-mass star forming regions such as 
the Taurus clouds, the typical stellar mass is of order 
0.5 Mq. As a rough approximation, we will limit the 
masses of all formed stars to the same value, M*. As 
soon as the stellar mass reaches M,, the accretion is ter- 
minated suddenly. At the same time, we eject the dense 
material surrounding the formed star in an outflow, with 
a speed set by the reqiiiremcnt that the outflowing gas 
carry away a certain amount of (linear) momentum (see 
also Allen & Shu 2000). We envision the outflowing gas 
as the core material swept up by a fast protostellar wind 
emanating from close to the central star in a momentum- 
conserving fashion. The wind momentum is given by 

= MwindKvind = 100 / {Mq km/s), (1) 

where the dimensionless parameter / is the product of 
the wind speed V^ind in units of 100 km/s and the ratio 
of the mass ejected in the wind M^ind to the stellar mass 
. The protostellar outflows included in our simulation 
provide a potential means to replenish the supersonic 
turbulence, which decays quickly despite the presence 
of a dynamically important magnetic fleld (Stone et al. 
1998; Mac Low et al. 1998; Padoan & Nordlund 1999). 
They may also directly induce star formation. 

As in Paper I, we will adopt a sheet-like cloud geom- 
etry, where ambipolar diffusion can be computed accu- 
rately and expeditiously. Our goal is to investigate the 
interplay between magnetic flelds and supersonic turbu- 
lence, taking into account the effects of the stars formed 
in the c;lou(l and their protostellar outflows. Specifically, 
we want to address two fundamental issues: the efficiency 
of star formation and the formation of quiescent cores 
out of a turbulent background medium. The remainder 
of the paper is organized as follows. In § 2, we describe 
our mathematical formulation of the problem. Numeri- 
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cal results are described and interpreted in § 3 through 
§ 5, and their implications are discussed in § 6. The last 
section, § 7, contains a brief summary. 

2. FORMULATION OF THE PROBLEM 
2.1. Governing Equations 

We consider clouds threaded by well ordered magnetic 
fields, and adopt the standard thin-sheet approximation, 
taking advantage of the tendency for strongly magne- 
tized clouds to settle along the field lines into a fiat- 
tened configuration in a dynamical time (e.g., Fiedler 
& Mouschovias 1993; Ostriker, Gammie & Stone 1999). 
We restrict the gas motions to the plane of matter dis- 
tribution, taken to be the x-y plane of a Cartesian sys- 
tem {x,y,z). Possible coupling between the motions in 
the plane and in the vertical direction will be addressed 
in future 3D simulations (see Krasnopolsky & Gammie 
2005). The 2D cloud evolution is governed by a set of 
vertically integrated MHD equations. These equations 
are the same as those used in our previous investigations 
(Li & Nakamura 2002; Nakamura & Li 2003; Paper I). 
They are described below for reference. 

The vertically integrated equation for mass conserva- 
tion is 

— -1-V-(EV)=.0, (2) 

where V = (Vx.Vy) is the velocity of the bulk, neutral 
cloud material in the plane of mass distribution. The 
momentum equation is 

where B =^ {B^, By) and B^ are the horizontal and verti- 
cal components of the magnetic field, and g = [g^, gy) is 
the horizontal component of the gravitational accelera- 
tion. The quantity P is the vertically integrated pressure, 
and H the half-thickness of the sheet. The latter is re- 
lated to the column and mass densities S and p through 



H = ^/{2p). 



(4) 



The equation governing the evolution of the vertical mag- 
netic field component is 



dB, 
dt 



V ■ (S,Vb) = , 



(5) 



where = (Vb,^, Ve.y) is the velocity vector of mag- 
netic field lines in the plane. 

The governing equations lO, and lO are supple- 
mented by the usual isothermal equation of state, which 
in a vertically integrated form becomes 



(6) 



where Cs is the isothermal sound speed. The isothermal 
condition ((HJ breaks down in low density regions where 
heating by UV radiation is important and in high density 
regions where thermal radiation becomes trapped; the 
deviation from isothermality is ignored at low densities 
since the thermal pressure is dominated by the turbulent 
pressure, and at high densities because the thermal pres- 
sure plays a minor role in the densest regions that are 
collapsing dynamically. 



The velocity of magnetic field lines in the induction 
equation (|3J| can be computed from 



Vb-v4 



B,B 
27r 



' R2 
47r 



(7) 



where t^ is the coupling time between the magnetic field 
and neutral matter. Realistic computations of are 
complicated by dust grains, whose size distributions are 
uncertain in dense clouds (e.g., Nishi, Nakano & Ume- 
bayashi 1991). For our purposes, we adopt the expression 



tr = 



1.4 



(8) 



which is valid in the simplest case where the coupling 
is provided by ions that are well tied to the field lines 
and the ion density pi oc p^/'^. For the coefficient C, we 
adopt the fiducial value 1.05 x 10"^ cm^/^g-^/^s"^ (Shu 
1991). The factor 1.4 in the above expression accounts 
for the neglect of ion-helium collision, whose cross sec- 
tion is small compared to that of ion-hydrogen collision 
(Nakano 1984; Mouschovias & Morton 1991). In rela- 
tively diffuse regions where the visual extinction is less 
than a few magnitudes, photoionization becomes impor- 
tant, which can enhance the magnetic coupling (McKee 
1989). An enhanced coupling at densities below the av- 
erage value adopted in our calculations does not change 
the results qualitatively. 

The mass density that appears in equations Q and 
^ can be determined from the condition of approxi- 
mate force balance in the vertical direction (Fiedler & 
Mouschovias 1993), which yields 



P 




(9) 



The two terms in the brackets represent, respectively, the 
gravitational compression and magnetic squeezing of the 
cloud material. The quantity Pe is the ambient pressure 
that helps confine the sheet, especially in low column 
density regions where gravitational compression is rela- 
tively weak. 

2.2. Initial and Boundary Conditions 

The equations @-® governing the cloud evolution 
are solved numerically. For initial conditions, we choose 
a uniform distribution of mass in the x-y plane, threaded 
by a magnetic field of constant strength Bq in the z direc- 
tion. In the absence of magnetic field, the cloud would be 
gravitationally unstable to perturbations of wavelength 
greater than the Jeans length Lj = Cg/(GSo) (Lar- 
son 1985). Here, Eq is the column density normal to 
the sheet. From the Jeans length Lj and sound speed 
Cs, we can define a time scale for gravitational collapse 



g _ Lj/cs = Cs/(GSo) (e.g., Ostriker et al. 1999). The 



Jeans length Lj and collapse time tg provide two basic 
scales for our problem. We adopt a square computation 
box of length L 3> Lj on each side, with periodic con- 
ditions imposed at the boundaries. The region under 
consideration thus contains a large number (~ \L/LjY) 
of thermal Jeans masses. Individual regions of size com- 
parable to Lj are expected to collapse to form stars or 
stellar systems on a time scale comparable to tg, unless 
other agents intervene. 
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The gravitational instability is suppressed if the field 
strength Bq is greater than a critical value 27rG^/^I]o (so 
that the cloud is magnetically subcritical) in the ideal 
MHD hmit (Nakano & Nakamura 1978). In a typical 
molecular cloud that is lightly ionized, the instability 
can still grow, albeit on the longer time scale of am- 
bipolar diffusion (Langer 1978). We have previously fol- 
lowed the nonlinear evolution of gravitational instability 
in weakly ionized, magnetically subcritical clouds, and 
found that ambipolar diffusion-driven fragmentation can 
readily lead to binary and multiple star formation in 
quiescent cloud cores in the presence of small perturba- 
tions (Nakamura & Li 2002, 2003; Li & Nakamura 2002; 
see also Indebetouw & Zweibel 2000 and Basu & Ciolek 
2004). A strong turbulence is expected to have a larger 
effect on cloud fragmentation. Some aspects of star for- 
mation in turbulent magnetized clouds were examined 
in Paper L They are investigated in greater depth in the 
present paper. 

To mimic the turbulent motions observed in molecular 
clouds, we introduce into the initially uniform cloud a su- 
personic velocity field at the beginning of the simulation. 
Following the usual practice (e.g., Ostriker et al. 1999; 
Klessen et al. 2000), we prescribe the turbulent velocity 
field in Fourier space, with a power-spectrum v'^ oc fc~". 
Unless noted otherwise, we will choose a power index 
oi n — 3, which is compatible with the Larson's (1981) 
size-dispersion relation in two spatial dimensions (Gam- 
mie 2003, priv. comm.). Numerically, we consider only 
discrete values ki = 2iTi/L [i = —N/2, N/2, with 
N being the number of grid points in each direction) for 
the two wavenumber components, kx and ky. For each 
pair of kx and ky, we randomly select an amplitude from 
a Gaussian distribution consistent with the power spec- 
trum for the total wavenumber k — ^k'^ + fc^ of that 

pair and a phase between and 27r. The resulting field 
is then transformed into the physical space to obtain the 
distribution of a velocity component at each grid point. 
The distributions of 14 and Vy are generated indepen- 
dently, and the final velocity field is scaled so that the 
root mean square (rms) Mach number of the flow A4 has 
a prescribed value. The velocity field so specified tends to 
have a strong compressive component. For comparison, 
we also investigate cases with incompressible turbulent 
velocity fields, generated by imposing the condition that 
the wavenumber vector be perpendicular to the direction 
of wave propagation (Ostriker et al. 1999). 

2.3. Dimensional Scalings 

The actual numerical computations are carried out us- 
ing non-dimensional quantities. The dimensional units 

we adopt are Cs = 1.88 x W^T^q"^ cm s~^ for speed 
(where Tio is the temperature in units of 10 K), Eq = 
4.68 X 10~^Av g cm~^ for surface density (where Ay 
is the visual extinction vertically through the sheet for 
standard grain properties), and Bq — 7. 59 AyTo fiG 
for magnetic field strength (where Fg is the fiux-to- 
mass ratio in units of the critical value). The units for 
length, time, and mass are, respectively, the Jeans length 
Lj = Cg/(GEo) = 0.37 Tio/Ay pc, gravitational collapse 

time tg = Cs/(GEo) = 1.90 x lO'^Tl^'^/Ay yr, and Jeans 
mass Mj = EqL} = 3.02 T^o/Ay Mq. We normalize 



the external pressure Pe by the "gravitational pressure" 
27rGEQ, and set its dimensionless value to 0.025 for all 
models. As long as the external pressure is significantly 
lower than the "gravitational pressure" , its exact value 
does not affect the cloud evolution much. 

2.4. MHD Code, Outflow Implementation and 
Lagrangian Particles 

We carry out numerical calculations using an MHD 
code that treats self-consistently the evolution of mag- 
netic field in a sheet-like geometry, including ambipo- 
lar diffusion. The code is described in Li & Nakamura 
(2002), modified here to take into account the periodic 
boundary conditions imposed in this problem. It was 
based on the code originally developed in Nakamura & 
Hanawa (1997). Briefly, we solve the hydrodynamic part 
using Roe's TVD (Total Variation Diminishing) method 
given in Hirsch et al. (1990). The vertical component 
of magnetic fleld is evolved in a way similar to the col- 
umn density, since both of them satisfy the same form 
of continuity equation (see equations and [S]). We 
define a magnetic potential outside the sheet, which is 
used to determine the magnetic forces on the sheet, in 
a manner similar to the determination of gravitational 
forces from the gravitational potential. Both the mag- 
netic and gravitational potentials are solved using a con- 
volution method based on FFT, as described in Hock- 
ney et al. (1981). The computational timestep is chosen 
small enough to satisfy both the usual Courant condition 
and the condition for treating ambipolar diffusion stably 
in an explicit code such as ours. We have verified that in 
the ideal MHD limit, the numerical diffusion of magnetic 
flux in our code is negligible, with the flux-to-mass ratio 
conserved within the machine roundoff error. 

We flx the size of our simulation box to 10 times 
the Jeans length, i.e., L = 10 Lj — 3.7 Tio/Ay pc. 
There are, therefore, 10^ thermal Jeans masses (or 
302 T^q/Av Mq) in our (two dimensional) computational 
domain. A 512 x 512 grid is used, with a cell size of 
AL = 7.2 X 10"3 Tio/Av pc or about 1500 Tio/Ay AU. 
Obviously, the grid is too coarse to resolve the flow pat- 
tern close to a forming star. Nevertheless, once dynami- 
cal collapse is initiated in a dense region, the mass accu- 
mulated in the highest density cell should provide a fair 
estimate of the stellar mass. As mentioned in the intro- 
duction, the stellar mass is probably limited by proto- 
stellar outflows, although the details of the infall-outflow 
interaction remain uncertain. For simplicity, we imple- 
ment into the MHD code the following procedure for stel- 
lar mass determination: as soon as the column density 
in a cell reaches a threshold Smax, we replace it with 
a much lower floor value S,„in. The mass extracted, 
= (Sniax — Sniin)(Ai)^, is placcd in a Lagrangian 
particle, which represents a formed star. Our canonical 
choices are Emax = 460So and S,„in = 30So, which yield 
a stellar mass Af* =0.5 Tiq/Av Mq. We keep the mag- 
netic flux in the cell unchanged, to reflect the fact that 
stars are born with a flux-to-mass ratio much smaller 
than the critical value. 

At the same time that the stellar mass is extracted 
from a cell, we set the mass Affiow within a region of 
radius 5AL = 3.6 x 10^^ Tio/Ay pc (centered on the cell) 
into a radial motion (in the plane of mass distribution) , 
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with a speed 

Vaow = A/indA^flow = 100 / 



km/s, (10) 



where equation (^) for protostellar wind momentum has 
been used. The outflow parameter / is the product of 
the wind speed in units of 100 km/s and the ratio of the 
mass ejected in the wind to the steUar mass. For opti- 
caUy revealed classical T Tauri stars, the wind speed is 
typically a few hundred km/s, and ~ 10% of the mass 
accreted through the disk is ejected in the wind (Cal- 
vet 1997). If the wind is launched magnetocentrifugally 
from a circumstellar disk (e.g., Konigl & Pudritz 2000), 
the wind speed is expected to be somewhat lower during 
the embedded, protostellar accretion phase than during 
the revealed phase, because of a lower stellar mass (and 
thus slower disk rotation) and a higher mass loading. If 
protostellar winds carry away a similar fraction of the 
accreted mass as the T Tauri winds, then 0.1, for a 
reasonable wind speed of 100 km/s. This is the canoni- 
cal value that we will adopt. A complication is that most 
of the wind momentum may escape in a direction per- 
pendicular to the mass distribution in 3D, which would 
lower the value of /. To allow for this possibility, we also 
consider cases with a much reduced / = 0.01. 

In our model, the termination of mass accretion and 
turn-on of outflow are assumed to be sudden. These can 
be justified on the ground that the bulk of mass accretion 
and outflow occur during the Class phase of low-mass 
star formation, which lasts for only a few times 10^ years 
(Andre et al. 2000), much less than the collapse time 

tg = 1.90 X yr at the average cloud density. 

Once created, the (stellar) Lagrangian particles are 
evolved in the collective gravitational potential of both 
the cloud material and the stars. To avoid singularities, 
we soften the gravitational potential of each star within 
a radius 3AL = 2.2 x 10"^ Tio/^v pc of the star. Al- 
though the choice of the radius is somewhat arbitrary, 
the gravitational softening may mimic the effects of a 
stellar wind, which may prevent a formed star from ac- 
creting additional ambient material as it moves around 
in the potential well of the cloud. 

3. EFFICIENCY OF STAR FORMATION IN TURBULENT 
MAGNETIZED CLOUDS 

3.1. Standard Simulation 

A number of parameters are needed to specify the tur- 
bulence, magnetic field and protostellar outflow in our 
problem. Some of these will be explored in subsequent 
subsections. Here, we focus on a standard simulation 
with a set of standard parameters to illustrate the basic 
results. It will also serve as a standard against which 
other simulations are compared. Since one of our em- 
phases is on ambipolar diffusion, which is expected to 
play a more important role in magnetically subcritical 
clouds than in supercritical clouds, we choose a "stan- 
dard" dimcnsionless flux-to-mass ratio Fq = 1.2. The 
cloud is thus moderately subcritical. At the beginning 
of computation t = 0, we stir the cloud with a random 
velocity field of power index n ~ 3 (so that the tur- 
bulent energy is dominated by large-scale motions) and 
rms Mach number Ai — 10, generated using the pro- 
cedure outlined in § 12.21 Unless noted otherwise, this 



"standard" random realization of turbulent velocity field 
is applied to other simulations. The star created in the 
collapse of a dense core is handled by a Lagrangian par- 
ticle according to the prescription outlined in § 12.41 The 
remaining core material is ejected in an outflow, with the 
outflow parameter / set to the canonical value 0.1. The 
parameters of the standard simulation and its variants 
are listed in Table 1. 

We have followed the cloud evolution in the standard 
simulation (Model SI in Table 1) well beyond the grav- 

1/2 

itational collapse time tg = 1.90 T^^ jAy Myr of the 
initial cloud. A movie of the evolution can be obtained 
from the authors. The main results are shown in Figs. 1- 
3. In Fig. 1, we plot the most important quantity that 
we are after, the star formation efficiency (SFE here- 
after), as a function of time. The first star forms around 

t = 0.65 tg = 1.24 tII'^/Av Myr. Thereafter, the num- 
ber of stars increases slowly but steadily. By the time 

i = 4 = 7.60 tII'^/Av Myr, 32 stars of O.hT^^/Av Mq 
each have formed, yielding an accumulative SFE of 5.3%. 

The low efficiency of star formation is the hallmark of 
the standard simulation. The efficiency is low because of 
the strong magnetic field imposed on the cloud initially. 
The field strength is such that the majority of the cloud 
mass remains magnetically subcritical even in the pres- 
ence of a supersonic turbulence. Pockets of supercritical 
material that are capable of forming stars are created 
through ambipolar diffusion, which is a relatively slow 
process. The ambipolar diffusion is governed by equa- 
tion IJj) , which can be cast into a more transparent form 



Vb-V 



0.17 



p2 



B 



where the auxiliary quantity 

^\Bl 



D = l + 



B'i 



(11) 



(12) 



is the column density enhancement factor due to mag- 
netic compression in the vertical direction (with the usu- 
ally small external pressure term in equation [H] ignored) . 
For a cloud near magnetic criticality F ~ 1, the drift 
speed of magnetic field lines relative to the bulk neutral 
material, Vb — V , is much smaller than the sound speed 
Cs in regions of mild field variation (where the field lines 
are nearly vertical, i.e., B^ + By <^ B^, and Bz varies 
gradually in space). The diffusion is faster in shocks 
where Bz changes over a shorter distance and the field 
lines are bent more strongly away from the vertical di- 
rection. The enhanced diffusion in shocks leads to accel- 
erated production of supercritical material in localized 
regions of the turbulent, magnetically subcritical clouds. 

The magnetically supercritical regions created through 
ambipolar diffusion are shown in Fig. 2. The figure dis- 
plays the snapshots of the standard cloud at six repre- 
sentative times t = 0.5, 1.0, 1.5, 2.0, 3.0, and 4.0 tg, with 
white contours separating the "barren" subcritical mate- 
rial that cannot form stars from the "fertile" supercrit- 
ical material that has the potential for star formation. 
At early times, supercritical material is created in nu- 
merous small pockets, as a result of the power spectrum 
of the initially imposed turbulent velocity field. The ve- 
locity distribution is such that regions of small character- 
istic lengthscale collide first, creating a large number of 
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shocklets (see also Gammie et al. 2003), where ambipo- 
lar diffusion is accelerated. The amount of mass involved 
in each shocklet is small, less than the "effective" Jeans 
mass, which is increased over the purely thermal value 
by both the magnetic pressure and partial cancellation 
of the gravitational force by the magnetic tension force 
(Nakamura & Hanawa 1997; Shu & Li 1997). These su- 
percritical regions do not collapse promptly. They are 
swept up at later times by the larger scale, converging 
turbulent motions into a network of more massive fila- 
ments. The supercritical filaments tend to be denser than 
the subcritical background because their self-gravity is 
stronger than the magnetic tension force, especially at 
late times. The cancellation of the tension force by self- 
gravity "softens" the supercritical regions, making them 
easier to compress by turbulence and harder to rebound 
after compression. 

The mass fraction of the supercritical material is plot- 
ted in Fig. 3 as a function of time. It increases sharply 
initially, as a result of the large gradients in the initial 
turbulent velocity field, which create strong shocklets 
that accelerate ambipolar diffusion. The second increase 
occurs around one turbulence crossing time, defined here 
as the time for the large-scale turbulence to cross half of 
the simulation box, 



For the standard simulation, the crossing time = 
0.5 tg. The increase is probably caused by the large-scale 
shock compression. Thereafter, the mass fraction stays 
more or less constant (roughly a third) over a period of 
several (average) gravitational collapse times. This value 
provides an upper limit to the efficiency of star forma- 
tion. 

Most of the supercritical filaments remain magnetically 
supported to a large extent. They last for a time long 
compared with the local dynamical collapse time. The 
longevity is unique to the shock-produced filaments that 
occur in magnetically subcritical clouds through ambipo- 
lar diffusion. Similar filaments created in weakly mag- 
netized, supercritical clouds are short-lived; they frag- 
ment and collapse promptly in a dynamical time. On 
the other hand, in the absence of ambipolar diffusion 
(i.e., ideal MHD limit), converging turbulent flows can 
still produce dense filaments in subcritical clouds. The 
filaments remain subcritical, however. They quickly dis- 
perse away once the external confinement disappears 
(Vazquez-Semedani et al. 2005). Evidently, the ambipo- 
lar diffusion in our standard simulation has weakened 
the magnetic field inside the filaments enough for them 
to "stick" together after compression, but not so much 
that their self-gravity would completely overwhelm the 
magnetic forces and collapse the filaments dynamically. 
The magnetic support of supercritical filaments, which 
contain a minor fraction of the cloud mass to begin with, 
further reduces the rate and efficiency of star formation. 

Stars are formed in dense cores of supercritical fila- 
ments. These cores will be discussed in depth in the 
next section, § 4. Here we merely point out that only 
a small fraction of the supercritical material resides in 
the densest regions that are most intimately linked to 
star formation. Fig. 3 shows that the mass fraction is 
about 10% at column densities above 10 times the aver- 



age value, typical of the low-mass dense cores identified 
in nearby dark clouds. The dearth of dense gas further 
limits the star formation efficiency. 

As we will see in § 4, some of the dense cores expand 
and disperse away without star formation. For those that 
do collapse to form stars, the local efficiency is limited 
by outflows. In our simulation, mass accretion from the 
core is cut off once the stellar mass reaches a pre-specified 
value (typical of the low-mass stars formed in nearby 
dark clouds) . The remaining core material is blown away 
in an outflow. Some of the cavities evacuated by outflows 
are evident in Fig. 2. On a global scale, the outflows can 
in principle stimulate further star formation by compress- 
ing supercritical regions into collapse, and by speeding up 
ambipolar diffusion in the material that they sweep up. 
They do not, however, fundamentally change the ineffi- 
cient nature of star formation in magnetically subcritical 
clouds such as the one in our standard simulation. 

3.2. Varying Outflow Strength of Standard Simulation 

The amount of protostellar outflow momentum that is 
coupled to the ambient medium is uncertain. In particu- 
lar, there is the possibility that most of the outflow mo- 
mentum escapes perpendicular to the plane of mass dis- 
tribution. To explore this possibility, we rerun the stan- 
dard simulation with the outflow parameter / reduced 
by a factor of 10, to 0.01 (Model S2 of Table 1). The 
results are shown in Fig. 4, where the snapshots of the 
cloud at six representative times t = 0.5, 1.0, 1.5, 2.0, 3.0 
and 4.0 tg are displayed, as in Fig. 2 for the standard 
simulation. Comparison of the two flgures shows that 
the cloud morphologies are rather similar, especially at 
early times, despite the factor-of-tcn difference in out- 
flow strength. The main difference lies in outflow cavi- 
ties which, as expected, are less prominent in the weaker 
outflow case. 

One may naively expect the stronger outflows in the 
/ = 0.1 case to stimulate more star formation. This turns 
out not to be the case. The number of stars formed near 
the end of the simulation a.t t — i tg is lower in the 
/ = 0.1 case (32) than in the / = 0.01 case (40) (see 
also Fig. 6 below). The reason, we believe, is that a 
stronger outflow tends to disperse the supercritical ma- 
terial around a formed star over a larger region, making 
it more difficult to re-condense. In any case, the differ- 
ence in the efficiency of star formation is modest (only 
25%), which leads us to believe that the SFE is relatively 
insensitive to the strength of protostellar outflow. 

3.3. Varying Turbulent Velocity Field of Standard 
Simulation 

The kinetic energy of the turbulence in molecular 
clouds is thought to be dominated by large-scale motions 
(Larson 1981). The exact velocity field is unknown, how- 
ever. It is usually modeled through random realization of 
a power spectrum specified in Fourier space, as described 
in § 12.21 To check the dependence of our results on the 
prescribed turbulent velocity field, we carried out a set 
of 6 simulations, keeping the rms Mach number fixed at 
= 10. The model parameters are listed in Table 1. 
Simulations Rl and R2 arc the same as the standard sim- 
ulation except that the turbulent velocity field of power 
index n = 3 is realized in a different way. Simulations II 
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and 12 also have power spectra of n = 3 but the veloc- 
ity fields are chosen to be incompressible. Simulations 
Nl and N2 have random velocity fields of power spectra 
with n = 5 and 1, respectively. The results are shown in 
Fig. 5. 

All four different random realizations of turbulent ve- 
locity field of the same n — 3 power spectrum yield a 
SFE close to that of the standard simulation. The SFE 
of the 71 = 5 case, where the power is more dominated 
by the large scales, is hardly distinguishable from those 
of the n = 3 cases. The n = 1 case, on the other hand, 
is quite distinct from all other cases. Its star formation 
did not start until about t = 3.5 tg, after most of the 
initial turbulence has already decayed away. The lower 
efficiency in this case is probably due to a larger fraction 
of the turbulent energy residing on small scales, which 
dissipates quickly without compressing enough mass to 
enable gravitational collapse. Based on this set of sim- 
ulations, we conclude that the low SFE we obtained in 
the standard simulation is robust. It is, we believe, a 
generic feature of star formation in turbulent, moder- 
ately subcritical clouds. We now turn our attention to 
non-magnetic and marginally magnetically supercritical 
clouds, which behave differently. 

3.4. Marginally Magnetically Supercritical and 
Non-magnetic Clouds 

Weaker magnetic fields lead to higher rates of star for- 
mation. We demonstrate this trend with two examples. 
Both have the same initial mass distribution and tur- 
bulent velocity field as in the standard simulation. The 
only difference is that their fiux-to-mass ratio Fq — 0.8 
(Model Ul) and (Model 01), respectively, instead of 
1.2. The results are shown in Figs. 6 and 7. We consider 
these two cases in turn. 

The marginally magnetically supercritical cloud (Fq = 
0.8) forms stars continuously over a time much longer 
than the turbulence crossing time {t^ = 0.5 tg), with a 
rate that decreases with time (Fig. 6). The decrease is 
due, at least in part, to the fact that, as more (weakly 
magnetized) stars form, the remaining gas becomes more 
magnetized. The average flux-to-mass ratio eventually 
exceeds the critical value and the cloud becomes subcrit- 
ical as a whole. Another reason may be the decay of tur- 
bulence, which is at a relatively low level despite the mo- 
mentum input from protostellar outflows (see §|Slbelow). 
Near the end of the simulation at t = 3.5 tg, 162 stars are 
produced, corresponding to a SFE of about 27%. The 
efficiency is in the range inferred for nearby embedded 
clusters (from - 10 to - 30%; Lada & Lada 2003). This 
similarity in SFE indicates to us that cluster-forming re- 
gions, although supercritical, may still be substantially 
magnetized. For comparison, we plot in the same figure 
the SFE for the non-magnetic cloud (Fq = 0). In this 
case, stars are formed much more quickly. More than 
half of the cloud material is converted into stars by the 
end of the simulation a,t t — 1.5 tg. The high efficiency is 
obtained despite the presence of a vigorous outflow from 
each of the formed stars. It is formally above the values 
given in Lada & Lada. 

As the outflow parameter decreases from the standard 
value / = 0.1 to 0.01, the efficiency of star formation 
becomes somewhat higher. Fig. 6 shows that, in the 
marginally supercritical (F = 0.8) cloud, about 10% more 



stars are formed by i = 3 in the / — 0.01 case than 
in the / = 0.1 case. The fractional increase is larger 
in the non-magnetic (Fq = 0) cloud. Its SFE in the 
weaker outflow case is even more inconsistent with the 
observationally inferred values (Lada & Lada 2003). 

Besides the rate and efficiency of star formation, the 
non-magnetic and marginally supercritical clouds differ 
in another important aspect: cloud mass distribution. 
The difference is illustrated in Fig. 7, where the snap- 
shots of the non-magnetic case at i = 0.5, 1.0, and 1.5 tg 
are shown in the first three panels and those of the mag- 
netized case ai t = 1.0, 2.0, and 3.0 tg in the last three. 
The contrast in cloud appearance is striking. Most of 
the gas in the non-magnetic cloud is concentrated in 
very thin dense filaments, which are clothed by little 
moderately dense envelopes. These essentially bare fil- 
aments are surrounded by large voids, similar to those 
commonly observed in pressureless cosmological simula- 
tions. In contrast, the filaments in the magnetized cloud 
are embedded within moderately dense gas, particularly 
at late times. A third difference is that the ffow speed 
is substantially higher in the non-magnetic cloud than in 
the magnetized cloud. This difference will be discussed 
further in § in connection with turbulence decay. 

3.5. Effects of Turbulent Mach Number 

The level of turbulence varies on different scales of 
molecular clouds and from place to place. Regions of 
stronger turbulence are expected to produce stars more 
vigorously. This trend is demonstrated in Fig. 8, where 
the SFEs of the standard simulation (with a turbulent 
Mach number — 10) and its two variants (Models S3 
and S4 in Table 1) of weaker turbulence [AA — 3 and 1, 
respectively) are plotted. As the initial turbulent Mach 
number decreases, star formation starts at a later time, 
with a longer quiescent period of little or no star forma- 
tion. In the case of A4 = 3, the first star forms around 

t = 1.3 tg (or - 2.5TII'^/Av Myr). It is followed by 
another non-star forming period of similar duration. Af- 
ter about t = 2.7 tg, the rate of star formation becomes 
comparable to, or even higher than the initially more 
turbulent standard case. The late increase in star for- 
mation activities is driven, we believe, by "spontaneous" 
ambipolar diffusion (as opposed to the "forced" ambipo- 
lar diffusion in strong shocks), which operates at all 
times, particularly in higher density regions. The "spon- 
taneous" ambipolar diffusion increases the mass fraction 
of supercritical material gradually after the initial strong 
turbulent compression. Even in this gentler mode of am- 
bipolar diffusion the turbulence plays a crucial role, by 
producing the density inhomogeneities needed for effi- 
cient diffusion. Once stars form, their outflows can fur- 
ther enhance the SFE by compressing nearby supercrit- 
ical regions (formed through ambipolar diffusion) into 
collapse. In this interpretation, the longest quiescent pe- 
riod, observed in the least turbulent Ai — 1 case, is 
related to the fact that its density peaks have the lowest 
values. Comparison of these three cases clearly demon- 
strates that in magnetically subcritical clouds star for- 
mation is accelerated by turbulence. 

Star formation in magnetically critical clouds is accel- 
erated by turbulence as well. The speedup of star forma- 
tion with increasing turbulent Mach number is shown in 
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Fig. 8, where the SFEs of a cloud identical to that in the 
standard simulation except for To — 1.0 (instead of 1.2) 
are plotted for three values of Mach number = 10, 3 
and 1 (corresponding to Models CI, C2, and C3 in Ta- 
ble 1, respectively). Note that, for a given M, the rate of 
star formation is higher for the critical (Fq = 1.0) cloud 
than for the standard (subcritical. Fq = 1.2) cloud. This 
trend is a clear demonstration that the star formation in 
a more strongly magnetized cloud is retarded to a larger 
extent. The behavior of star formation in critical clouds 
is intermediate between those of subcritical and super- 
critical clouds discussed in the previous subsections. 

4. DENSE CORES IN STANDARD SIMULATION 

Dense cores are the basic units of low-mass star for- 
mation. The best observed cores are those selected from 
nearby dark clouds of dispersed star formation that arc 
large enough for detailed mapping. Individual cores in 
cluster-forming regions such as p Oph tend to be smaller 
and thus less well resolved (Motte, Andre & Neri 1998). 
In this section, we will concentrate on the cores formed 
in our standard simulation, which we believe is represen- 
tative of the clouds undergoing dispersed star formation 
of low efficiency. 

In our standard simulation, dense cores are formed in 
supercritical filaments which, in turn, are produced in 
strongly shocked regions where ambipolar diffusion is ac- 
celerated. Since large-scale converging motions are can- 
celed out in post-shock regions, we expect the cores to de- 
velop in a relatively quiescent environment (e.g., Padoan 
et al. 2001). On the other hand, the cores and the fil- 
aments in which they are embedded are pushed around 
constantly by the ambient turbulent flow and occasion- 
ally by protostellar outflows, so perfectly quiescent cores 
are not to be expected. Furthermore, inside individual 
cores, self-gravity can in principle generate substantial 
infall motions. The velocity fields of some starless cores 
have been mapped and more will become available in the 
ALMA era. They provide a contact between theory and 
observations. 

There is no unique definition of dense cores. Ob- 

servationally derived properties of cores depend on the 
molecules or dust emission used to trace them. We adopt 
the following procedure to identify a core. First, we lo- 
cate all local maxima of the column density distribution 
above a threshold Ec = IOSq, where Eq is the average 
column density. Around each maximum, we then draw a 
set of density contours and identify the region enclosed 
within a level Se = 6 Sq as a dense core. The boundary 
value Ee is chosen to be roughly half of the peak den- 
sities of those cores that are barely above the threshold 
Ec. It corresponds to a visual extinction of magnitude 
6 (for the canonical choice of average Ay = 1), typical 
of many cores of dark clouds selected from the Palomar 
Sky Atlas (Myers et al. 1983). 

Cores are created, distorted and destroyed continu- 
ously. For illustration, we will concentrate on those iden- 
tified a,t t — 2 tg, the middle point of our standard sim- 
ulation. By this time, the initial turbulent velocity field 
has more or less relaxed. We identified 10 cores in total 
(excluding a partial ring created by the outflow of a newly 
formed star). There are double-cores sharing a common 
envelope; they are split into two single cores. The cores 
are plotted and numbered in Fig. 9, and shown in greater 



detail in Fig. 10. Plotted in each panel of Fig. 10 are 
the (column) density contours and vectors of the veloc- 
ity field relative to the density peak. The relative velocity 
is a measure of the (nonthermal) velocity dispersion of 
the core. Some properties of the cores are summarized 
in Table 2. In this table, the core radius (column 3) 
is evaluated as the radius of a circle enclosing the same 
area as the E^ = 6 Sq contour; for double cores, the area 
is split evenly between the two components. The aver- 
age flux-to-mass ratio in column 6 is defined as the total 
magnetic flux divided by the core mass, and the velocity 
dispersion in column 7 as the mass-weighted rms value 
of the flow velocity relative to the density peak. 

Prom Fig. 10 and Table 2, we find that the majority of 
cores have a subsonic internal velocity dispersion. One 
exception is core No. 2 (panel b in Fig. 10). Its peak den- 
sity of 118 Eq is the highest of all cores. The core has al- 
ready collapsed, with a "stellar" mass of 0.13 T^q/Av Mq 
located within the highest density cell. In effect, it con- 
tains a "Class 0" protostar still actively accreting mass 
from the surrounding massive envelope (Andre et al. 
2000); the large velocity dispersion of the core is thus 
gravitational in origin. Its mass accretion is terminated 

in another ~ 7 x IQ'^tH'^ /Ay yrs, comparable to the du- 
ration of the Class phase. The core that most closely 
resembles the collapsed core is core No. 6 (panel g). It 
has the second highest peak density of 41.2 Eq. System- 
atic infall has already developed in this core, although 
still at a subsonic speed. It is well on its way to star 
formation, formally producing a "finished" star in about 

0.18 tg = 0.34 tII'^/Av Myrs. The infaU motions are in 
contrast to the systematic outward motions in core No.l 
(panel a). This core has a peak density (11.3 Eq) slightly 
above the threshold Ec = 10 Eq. It is in the process of 
dispersing away. There are three other dispersing cores. 
They are cores No. 4, 5 and 10 (panels d, e and 1, respec- 
tively). These cores were recently formed out of material 
swept up by protostellar outflows, and tend to have rela- 
tively high internal velocity dispersions. The dispersing 
cores also have higher-than-average flux-to-mass ratios, 
with values between 0.81 and 1.03. The high level of tur- 
bulence and high degree of magnetization are probably 
the reasons why these cores are dispersing rather than 
condensing. 

The core with the lowest degree of magnetization and 
the lowest velocity dispersion is core No. 9 (panel k). This 
roundish core has the smallest mass and radius of all 
cores, and is apparently close to a mechanical equilib- 
rium. It lasts for nearly 0.45 tg (or about 5 times the 
local collapse time), before being induced to collapse by 
an protostellar outflow. 

The cores No. 7 (panel i) and 8 (panel j) are enclosed 
within a common envelope. Both have relatively high 
flux-to-mass ratios (0.91 and 0.80, respectively), which 
is probably the reason why they lived for a relatively 
long time (compared to their local collapse times) be- 
fore merging and forming stars around t = 2.36 tg. The 
longest living core is core No. 3. It can be traced to the 
end of simulation, for more than 20 times its local col- 
lapse time. At the time t = 2 tg (shown in panel c). the 
core has a relatively large velocity dispersion (close to 
the sound speed) and is being compressed by a converg- 
ing flow. The compression was not able to push the core 
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directly into collapse, probably because of its relatively 
large flux-to-mass ratio (0.88). This elongated core ap- 
pears to be initially far out of equilibrium. It gradually 
relaxes toward a round configuration. At later times, a 
rapid rotation develops. Rotational support is probably 
the main reason for its extraordinary longevity. 

Panels (f)-(h) of Fig. 10 illustrate the evolution of the 
collapsing core No. 6. As is typical, this core is initially 
elongated, as can be seen in panel (f) (at t = 1.78 tg), 
when it is assembled by a converging flow. At this early 
stage of formation, the core (as outlined by the density 
contour at Sg = 6 Sq) is extended and turbulent, with 
a supersonic velocity dispersion. It becomes rounder 
and more quiescent at later times, with a velocity dis- 
persion reaching down to 30% of the sound speed at 
t = 2 tg (shown in panel g). The slow, subsonic con- 
traction observed at this time soon gives way to a more 
rapid, supersonic collapse, which is shown in panel (h) 
(at t = 2.06 tg). The dip in velocity dispersion appears 
to be a common feature of star-forming cores during the 
period after the external compression is relaxed but be- 
fore dynamical collapse sets in. From panel (f) to (g), 
the core radius decreases by about a factor of two. The 
fact that this core is more compact at later times is con- 
sistent with the H^'^CO observations of Mizuno et al. 
(1994), which show a tendency for more evolved cores 
to be smaller in size. Note also that the density peak is 
offset from the geometric center of the core at all three 
times. Such a cometary shape is often observed in dense 
cores (Crapsi et al. 2005). 

In summary, of the 10 cores shown in Fig. 9, one formed 
a star already (No. 2), another is well on its way to star 
formation (No. 6). After being strongly perturbed, three 
more cores collapse to form stars (No. 7, 8, and 9). Of 
the remaining 5 cores, four expanded and disappeared as 
their peak densities drop below the threshold for a core 
(No.l, 4, 5, and 9). The fate of the last core (No.3) is 
uncertain. 

5. TURBULENCE DISSIPATION 

Supersonic turbulence is known to dissipate quickly. 

This is also true in our simulations. The turbulence de- 
cay is illustrated in Fig. 11, where we plot the mass- 
weighted rms Mach number of the gas 
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of the standard simulation as a function of time. The 
Mach number drops from the initial value 10 to ~ 4 in 
one turbulence crossing time (fx = 0.5 tg). The decline 
slows down at later times, with a factor of two reduction 
from one to eight crossing times. The slowdown results 
from a decrease in shock strength, as the flow speed ap- 
proaches the fast magnetosonic speed of the cloud. As 
a result of magnetic cushion, moderately supersonic mo- 
tions can survive for a long time, especially in relatively 
low density regions. 

For comparison, we also plotted the rms Mach number 
of the stars in Fig. 11. Except for an initial dip, the 
rms stellar Mach number increases gradually from ~ 2 to 
^ 3.5, before dropping back to 3. The Mach number is 
somewhat higher than that of the gas at late times. The 
most likely reason for this difference, we believe, is that 



stars experience the full gravity of all the gas inside the 
simulation box, whereas the gas feels only a magnetically 
diluted gravity. 

Protostellar outflows appear to have little effect on the 
replenishment of turbulence in our standard simulation. 
When the outflow strength is reduced by a factor of 10, 
the rms Mach number of the gas is hardly changed (see 
solid lines in Fig. 12). The ineffectiveness stems from the 
slow rate of star formation in this moderately subcritical 
cloud, where a relatively small amount of momentum is 
injected into the cloud over a long period of time. Even 
in the case of marginally supercritical cloud, where the 
star formation rate is much higher, the dissipated turbu- 
lence is not replenished significantly, as shown in Fig. 12 
(dotted lines). Only in the non-magnetic (Tq = 0) and 
strong outfiow (/ = 0.1) case, the rate of star formation is 
high enough that momentum injection from outflows can 
roughly offset dissipation and maintain the turbulence at 
a level comparable to the initial value for a collapse time. 
The high level of turbulence is also evident in the velocity 
fleld plotted in the first three panels of Fig. 7. When the 
outflow parameter / is reduced to 0.01, the turbiflence in 
the non-magnetic cloud decays with little replenishment, 
despite an increase in the rate of star formation. 

6. DISCUSSION 

6.1. Dispersed Star Formation of Low Efp,eieney 

Magnetic regulation was originally envisioned for the 
dispersed mode of inefficient star formation (Shu et al. 
1987). The best studied example of this mode is the Tau- 
rus molecular cloud complex. With a mass ~ 2 x 10^ M0 
(Ungerechts & Thaddeus 1987) and ^ 200 low-mass 
young stellar objects (Kenyon & Hartmann 1995), the 
overall star formation efficiency is about 1%. Palla & 
Stabler (2002) analyzed the star formation in Taurus 
both in space and in time. They found that it started 
~ 10 Myrs ago at a low level, and the pace of star forma- 
tion accelerated in the last few million years. Since the 
clouds must be at least as old as the oldest stars, Palla 
& Stabler concluded that the Taurus clouds are prob- 
ably older than ~ 10 Myrs. If true, the clouds would 
have lived for at least 5 times the gravitational collapse 
time tg, since their average visual extinction Ay ~ 1 
(Arce & Goodman 1999) and temperature T ~ 10 K 
yield tg ^ 2 Myrs. 

Hartmann (2003) questioned the above estimate of 
Taurus cloud lifetime. He pointed out that the oldest 
stars in Palla & Stabler sample tend to be stars more 
massive than ^Mq and that the ages of such stars are 
affected by uncertainties in birthline location. If these 
stars are discarded, the ages of the oldest stars (and thus 
the lower limit to the cloud lifetime) would be closer to 
5 Myrs (or ~ 2.5tg) than to 10 Myrs. The shorter lifetime 
was taken as evidence against the traditional, quasi-static 
picture of magnetically regulated star formation. In the 
presence of a supersonic turbulence, significant star for- 
mation can occur in moderately subcritical clouds in a 
few collapse times, so that even the shorter cloiid lifetime 
is not a problem for the picture of star formation accel- 
erated by turbulence and regulated by magnetic fields 
through ambipolar diffusion (for example, the SFE of the 
standard model reaches about 1 % in 1-1.5 or 2-3 Myrs 
for Fo = 1.2 and M = 10). 

Hartmann (2003) suggested that star formation in Tau- 
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rus is rapid and dynamic. A specific scenario envisioned 
was molecular cloud formation and prompt star forma- 
tion in converging flows of atomic clouds (Hartmann et 
al. 2001). It remains to be demonstrated, however, that 
only one percent of the molecular gas can be turned into 
stars in this scenario. Numerical simulations of large- 
scale turbulent flows containing many Jeans masses in- 
variably show rapid conversion of gas into stars, with 
efflciencies much higher than one percent (e.g.. Bate et 
al. 2003; see also section §1^31 • The high efflciency prob- 
lem is one of the main factors that motivated the original 
proposition that star formation in regions like Taurus is 
magnetically regulated (Shu et al. 1987). 

The combination of a strong magnetic field and super- 
sonic turbulence allows for a range of possibilities. In our 
standard simulation, the cloud is moderately subcritical 
(with a dimensionless flux-to-mass ratio F = 1.2) and the 
turbulence is relatively strong (with a rms Mach number 
A4 — 10). By the end of 4 gravitational collapse times 

{t = Atg = 7.60 T^l'^/Av Myr), 5.3% of the cloud mass 
has been converted into stars. It corresponds to an aver- 
age rate of star formation per unit mass 7 x 10~^Av /T^q^ 
yr~^. Interestingly, this (specific) rate is not far from the 
Galactic average of ^ 3 x 10~^ yr~^ for a star formation 
rate of 3Mq per year inside the solar cycle (Scalo 1986) 
and a molecular mass of IO^Mq in that region (Blitz 
& Williams 1999). One can make the star formation 
more inefficient by increasing the magnetic field strength 
and/or reducing the turbulent Mach number. Reducing 
the Mach number of the standard simulation from 10 to 
3, for example, lowers the SFE a,t t = A tg from 5.3% 
to 2.0%. The reduction factor is higher at earlier times, 
when few stars are formed in the lower Mach number 
case. The pace of its star formation accelerates at later 
times (see Fig. 8), as inferred by Palla & Stabler (2002). 

To summarize, we believe that dispersed star forma- 
tion of low efficiency in regions like the Taurus molecular 
clouds naturally occurs in moderately magnetically sub- 
critical clouds, through ambipolar diffusion that is ac- 
celerated by turbulence. Whether the bulk of the cloud 
material is indeed magnetically supported remains to be 
determined through Zecman measurements. Crutcher & 
Troland (2000) measured the line-of-sight field strength 
toward the LI 544 region in Taurus. The inferred de- 
projected flux-to-mass ratio F ~ 0.5 is in agreement 
with core formation models involving ambipolar diffu- 
sion (Basu & Mouschovias 1994; Nakamura & Li 2003; 
see the ratios for collapsing cores in column 6 of Ta- 
ble 2). For the general background of the Taurus clouds 
where Ay ^ 1 (Arce & Goodman 1999), the critical field 
strength is ~ 7.6 /zG. The line-of-sight component is ex- 
pected to be smaller than this value by a typical factor 
of '--^ 2 — 4 (depending on the cloud geometry; Shu et 
al. 1999). Whether such a weak field can be directly 
measured remains to be seen. Wolleben & Reich (2004) 
inferred a field strength exceeding ~ 20/iG from model- 
ing polarization measurements at 21 and 18 cm toward 
the Taurus cloud complex. The value need confirmation 
from Zeeman Measurements. 

6.2. Quiescent Cores and Inefficient Star Formation 

Starless cores define the initial conditions for low-mass 
star formation (Andre et al. 2000). Their structure pro- 



vides clues to how they are formed. Dust continuum ob- 
servations have shown that these cores have flat-topped 
density profiles (Ward-Thompson et al. 1994). This fea- 
ture can, however, be explained by essentially any type of 
core formation scenario. More discriminating are molec- 
ular line observations, which are sensitive to the internal 
motions of cores. Lee et al. (1999) and Gregersen & 
Evans (2000) carried out line surveys of nearby starless 
cores. A general conclusion is that optically thick lines 
can be displayed either to the blue or red side of optically 
thin lines, with an overabundance of sources with blue 
asymmetry. The simplest interpretation is that cores 
with "blue" asymmetry are contracting, and those with 
"red" asymmetry are expanding. Both collapsing and 
dispersing cores are observed in our simulations. Our col- 
lapsing cores tend to have higher column densities than 
the dispersing cores (see Table 2 for examples). This 
trend is consistent with the finding of Lee et al. (1999; 
see also Gregersen & Evans 2000) that the best infall 
candidates tend to have the highest column densities. 

There are two important characteristics for the dense 
cores formed in our standard simulation, which is repre- 
sentative of inefficient, dispersed star formation in mod- 
erately magnetically subcritical clouds. First, most cores 
are supercritical. This is the case, for example, for 9 out 
of the 10 cores shown in Fig. |2| the remaining one is 
very close to critical, with an average flux-to-mass ratio 
F — 1.03. The dearth of subcritical cores is in agree- 
ment with Nakano (1998), who argued that such cores 
are unlikely to exist because they require a strong exter- 
nal confinement to stay in equilibrium. Even if external 
confinement is available, it would be difficult to maintain 
the observed level of nonthermal velocity dispersion for 
a significant fraction of the core's (long) lifetime. These 
difficulties have led Nakano (1998) to abandon the pic- 
ture of core formation through ambipolar diffusion in fa- 
vor of one involving turbulence decay. We believe the 
abandonment is premature. We have shown that con- 
verging flows in a strong turbulence can create and con- 
fine high density subcritical regions long enough for accel- 
erated ambipolar diffusion to produce supercritical ma- 
terial, particularly in moderately subcritical clouds where 
the magnetic forces that drive the expansion are already 
offset by self-gravity to a large extent (so that the confine- 
ment requirement is much reduced; Shu & Li 1997) and 
a relatively small reduction in magnetic flux can make a 
piece of cloud supercritical. Once created, supercritical 
cores can be held together by self-gravity and are able 
to collapse to form stars more quickly than externally- 
confined subcritical cores (so that excessive turbulence 
decay is less of a problem). In our picture, the bulk of 
magnetic flux reduction occurs dynamically (in shocks), 
rather than statically, through ambipolar diffusion when 
the cloud is strongly turbulent. As the level of turbulence 
decreases, quasi-static ambipolar diffusion becomes more 
significant. 

The second important characteristic of the cores is that 
they have predominantly subsonic internal motions (see 
Table 2 for examples). This is in agreement with the 
finding that the velocity dispersions of starless cores in 
regions of inefficient, dispersed star formation such as the 
Taurus clouds tend to be thermally dominated (Jijina et 
al. 1999; Lee et al. 1999). Quiescent cores are difficult to 
arrange in turbulent, non- or weakly-magnetized clouds 
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that contain many Jeans masses: these cores expand or 
collapse dynamically unless the core mass happens to 
be close to the local thermal Jeans mass. In our simula- 
tions, quiescent cores appear to form out of the turbulent, 
moderately subcritical background without fine tuning. 
The reason, we believe, is that our cores, although su- 
percritical in general, remain strongly magnetized, with 
a dimensionless flux-to-mass ratio F close to unity (the 
critical value; see Table 2). Roughly speaking, the mag- 
netic tension force is — times the self gravity (Shu & Li 
1997, their equation [2.16]; the ratio would be exact if F 
is spatially constant; sec also Nakano 1998). Inside each 
core, the self gravity is balanced out by the magnetic 
tension force to a large extent (the "effective" gravity 
is 1 — F^ times the undiluted value). Barring occasional 
strong external perturbations, the magnetically levitated 
cores tend to contract or expand in "slow motion" regard- 
less of the amount of mass involved. The slower motions 
enable the quiescent cores to live longer than their dy- 
namically evolving counterparts, which tend to increase 
their number at any given time. The long-lived cores 
that eventually form stars have masses greater than the 
local thermal Jeans mass but are gravitationally stable 
due to significant magnetic support. They tend to oscil- 
late with subsonic speeds before rapid contraction. The 
oscillation, we speculate, may excite the type of pulsating 
flow pattern observed in B68 (Lada et al. 2003). Such 
cores start to collapse either after merging with other su- 
percritical regions or after significant additional flux loss 
due to ambipolar diffusion. As the level of turbulence de- 
creases, the latter becomes increasingly more important. 

Magnetic Icvitation also lies at the heart of the low 
efficiency of star formation. In a cloud that is subcritical 
overall, the self gravity is completely offset by magnetic 
forces except in localized pockets — the supercritical (F < 
1) filaments (see Figs. 2 and 4), which contain a fraction 
of the cloud mass (about 1/3 in our standard simulation). 
Only a small fraction of the supercritical material resides 
in dense cores (about 1/10 in our standard simulation) 
that are intimately related to star formation. Not all 
cores form stars; some of them disperse away without 
star formation. For those cores that do form stars, their 
efficiencies arc further limited by protostellar outflows, 
which can unbind the majority of core material. 

Our picture of turbulence-accelerated, magnetically 
regulated star formation has implications on astrochem- 
istry. Cores formed directly out of the low-density back- 
ground from rapid turbulent compression may have had 
little time for chemical reactions. They may be identi- 
fied with chemically "young" cores such as L1521E that 
are rich in carbon-chain species like CCS and show little 
evidence of CO depletion (Hirota et al. 2002; Tafalla & 
Santiago 2004) . Some of the dense cores expand back to 
low densities, and it would be interesting to explore the 
effects of these dispersed cores (where chemical deple- 
tion, for example, may have occurred) on the chemistry 
of the more diffuse region (Garrod et al. 2005). Others 
evolve further towards star formation, through contin- 
ued ambipolar diffusion and/or external perturbations 
(such as merging and shock compression). They may 
be identified with more evolved cores such as L1544 and 
L1521F (e.g., Crapsi et al. 2004) that are thought to be 
on the brink of star formation. The interplay between 
turbulence and magnetic fields in the presence of am- 



bipolar diffusion allows for a variety of possibilities for 
the evolutionary history of starless cores, which may be 
constrained through astrochemical modeling. 

6.3. Implications on Cluster Formation 

The majority of stars are thought to form in clusters. 
One of the best-studied nearby cluster-forming regions is 
the p Oph core. It has an estimated stellar mass 53 Mq 
and core mass 550 Mq, yielding a SFE of ~ 9% (Lada & 
Lada 2003). This efficiency may increase with time, to 
a value of order ~ 30%, the maximum value inferred for 
nearby embedded clusters (Lada & Lada 2003). SFEs of 
only a few tens of percent are not easy to obtain in tur- 
bulent cores that are not magnetized, unless the cores are 
constantly driven on a small enough scale (Klessen et al. 
2000). In our simulations, the non-magnetic cloud col- 
lapses promptly, forming stars with an efficiency between 
~ 55% to 70% at 1.5 collapse times, despite the pres- 
ence of protostellar outflows. The efficiency is reduced 
substantially as the field strength increases toward the 
critical value. The marginally supercritical cloud in our 
models has a field strength 20% below the critical value. 
Its star formation efficiency approaches ~ 30%, close to 
the observationally inferred maximum value. When and 
how the star formation is actually stopped in a Taurus- 
like or p Oph- like cloud remains uncertain (e.g., Palla & 
Stabler 2002). 

The non-magnetic cloud that we modeled has an- 
other difficulty. It produces a mass distribution that is 
dominated by a network of dense, thin filaments, with 
large voids in between, similar to the structure obtained 
in pressureless cosmological simulations (see Fig. 7). 
Nearby cluster-forming regions such as p Oph (Motte et 
al. 1998) and Serpens (Olmi & Testi 2002) cores appear 
to have a smoother mass distribution, with less cloud 
mass concentrating in dense clumps than in inter-clump 
regions (Johnstone et al. 2004). This difficulty disap- 
pears in marginally supercritical clouds, where the dense 
cores are surrounded by extended, moderately dense en- 
velopes. The better agreements in SFE and cloud mass 
distribution in the presence of a dynamically-important 
magnetic field lead us to believe that magnetic fields 
play a role even in cluster formation in the present-day 
Galaxy. This contention is consistent with available Zee- 
man measurements, which show that the (deprojected) 
flux-to-mass ratios are not far from the critical value 
in highly turbulent regions of high-mass star formation 
(Crutcher 1999). 

6.4. Caveats and Future Work 

Our calculations are two dimensional. The adopted 

sheet-like geometry is valid for magnetically-supported, 
multi- Jeans mass clouds that are subcritical or 
marginally supercritical but not necessarily for weakly 
magnetized clouds, unless they are flatten externally (by 
shocks, for example). Vertical force balance is assumed 
at all times, which may break down in shocked regions 
immediately following strong compression. Such regions 
can puff up temporally, lowering the volume density (and 
thus the rate of ambipolar diffusion) . The reduction may 
lead to a somewhat lower rate of star formation. Ad- 
dressing this issue would require 3D calculations. 

3D calculations will also be needed to address the is- 
sue of magnetic braking. Some of the cores formed in 
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our simulations rotate rapidly. A particular example is 
the core No. 3 shown in Fig. 9. It appears to be rota- 
tionally supported at late times, lasting for more than 
20 times the local gravitational collapse time. The rapid 
rotation is expected to be reduced by magnetic brak- 
ing (Basu & Mouschovias 1994), which is not included in 
our simulation. The braking should reduce the rotational 
component of the velocity field inside a core in general, 
perhaps making the core more quiescent. The reduction 
factor will probably be less than that obtained by Basu 
& Mouschovias for quasi-statically evolving clouds, since 
in our picture the core formation time is shortened by 
external compression. 

In our simulations, stars arc taken to have the same 
final mass, which is obviously an oversimplification. An 
alternative is to assume that the time for a star to accrete 
its mass is fixed. Some justification for this assumption 
comes from the finding that the range in accretion time 
is much narrower than that in stellar mass (by an or- 
der of magnitude; Myers & Fuller 1993). Resolving the 
accretion flow around each forming star would require 
a spatial resolution higher than that of our simulations. 
Future higher resolution simulations may allow us to di- 
rectly tackle the problem of initial mass function. We 
have assumed for simplicity that the magnetic fltix-to- 
mass ratio is initially uniform throughout the cloud. This 
assumption will be relaxed in future studies. 

The replenishment of turbulence is a potential prob- 
lem. In our standard simulation of inefficient star for- 
mation in a magnetically subcritical cloud, the velocity 
dispersion drops quickly in the first turbulence crossing 
time. The decline rate slows down at later times, with the 
flow velocity approaching twice the sound speed for the 
cloud as a whole; it is somewhat higher in the lower den- 
sity regions, where large patches of more or less coherent 
velocity field exist. More organized velocity fields tend 



to be less dissipative, as suggested by Palla & Stabler 
(2002). The typical flow velocity of a few times the sound 
speed is not as large as the velocity dispersions observed 
on the scale of cloud complexes. How the large-scale tur- 
bulence is generated and maintained remains uncertain 
(see Elmegreen & Scalo 2004 for a review). 

7. SUMMARY 

We have explored numerically the interplay of mag- 
netic fields and turbulence in the presence of ambipo- 

lar diffusion in a sheet-like geometry, taking into ac- 
coimt protostellar outflows in an approximate manner. 
In strongly turbulent, moderately magnetically subcriti- 
cal clouds, we found a dispersed mode of star formation 
with efficiencies at several percent level. Dense cores 
formed in such clouds tend to be quiescent, with sub- 
sonic internal motions, in agreement with observations. 
For strongly turbulent, marginally supercritical clouds, 
the efficiency increases to a few tens of percent, which is 
in the range inferred for nearby cluster-forming regions. 
If not regulated by magnetic fields at all, the efficiency 
of star formation exceeds half in a collapse time, despite 
vigorous protostellar outflows. We conclude that a happy 
marriage between turbulence and magnetic fields must be 
sought, particularly for the inefficient, dispersed mode 
of star formation, perhaps for the more efficient, clus- 
tered mode as well. Further refinements of this picture of 
turbulence-accelerated, magnetically regulated star for- 
mation require 3D simulations. 
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TABLE 1 
Model Parameters 



Model 


To 


M 


n 


/ 


Turbulence " 


Notes 


SI 


1.2 


10 


3 


0.1 




.5 LdliLldl LI lilULiei, & UULl ILlLdl 


32 


1.2 


10 


3 


01 






Rl 


1.2 


10 


3 


1 




t_lllit:l eiili 1 CLlH_lUlli i ecLilZiCLLltJli 


R2 


1.2 


10 


3 


0.1 


Q 


ri 1 frV^rfint vctTirinTn i'fmli^ctl"ir\Ti 




1 9 


1 n 

-LU 


Q 

o 


n 1 


J) 


diflFerent rjmdom retilizfitiou 


12 


1.2 


10 


3 


0.1 


E 


difTerent random realization 


Nl 


1.2 


10 


5 


0.1 


F 


different power spectrum 


N2 


1.2 


10 


1 


0.1 


G 


different power spectrum 


Ul 


0.8 


10 


3 


0.1 


A 


marginally supercritical 


U2 


0.8 


10 


3 


0.01 


A 


supercritical, weak outflow 


Ol 


0.0 


10 


3 


0.1 


A 


non-magnetized 


02 


0.0 


10 


3 


0.01 


A 


non-magnetized, weak outflow 


S3 


1.2 


3 


3 


0.1 


A 


intermediate turbulence 


S4 


1.2 


1 


3 


0.1 


A 


weak turbulence 


CI 


1.0 


10 


3 


0.1 


A 


magnetically critical 


C2 


1.0 


3 


3 


0.1 


A 


critical, intermediate turbulence 


C3 


1.0 


1 


3 


0.1 


A 


critical, weak turbulence 



^The letter A denotes the random realization of turbulent velocity fleld used in the standard simulation, and B, C, D and E for different 
random realizations of the same power spectrum as A. The velocity flelds in D and E are incompressible. The power spectrum in E (F) is 
steeper (shallower) than that in A-E. 



TABLE 2 

Properties of the Cores Identified in the Standard Simulation at t = 2 tg " 



Core 


Core 


Core 


Peak 


F at 


Average 


Velocity 


Notes on Core 


No. 


Mass 


Radius 


S 


E Peak 


F 


Dispersion 


Evolution 


1 


0.32 


0.11 


11.3 


0.85 


0.89 


0.48 


expanding, dispersed 


2 


0.87 


0.15 


118 


0.60 


0.67 


1.40 


collapsing, star formed at t=2.04 tg 


3 


0.32 


0.10 


18.0 


0.79 


0.88 


0.97 


fast rotation develops, fate unknown 


4 


0.48 


0.12 


21.2 


0.82 


0.85 


1.41 


expanding, dispersed 


5 


0.35 


0.12 


12.9 


0.77 


0.81 


0.64 


expanding, dispersed 


6 


0.91 


0.17 


17.4 


0.85 


0.90 


1.33 


core No. 6 at t=1.78 tg 


6 


0.34 


0.095 


41.2 


0.57 


0.61 


0.30 


collapsing, star formed at t=2.18 tg 


6 


0.33 


0.077 


90.2 


0.53 


0.78 


1.18 


core No. 6 at t=2.06 tg 


7 


0.32 


0.13 


11.8 


0.85 


0.91 


0.55 


collapsing'', star formed t=2.36 tg 


8 


0.42 


0.13 


16.9 


0.68 


0.80 


0.40 


collapsing'^, star formed t=2.36 tg 


9 


0.11 


0.068 


11.4 


0.53 


0.57 


0.15 


collapsing'^, star formed t=2.44 tg 


10 


0.54 


0.15 


11.7 


0.85 


1.03 


0.80 


expanding, dispersed 



^'the units for the dimensionless quantities listed can be found in S 12.31 
^this core merges with core No. 8 to form a star, 
'^this core merges with core No. 7 to form a star. 

"^this core collapses to form a star after colliding with a dense region located between cores No. 7 and 9. 
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Fig. 2. — Snapshots of the standard simulation at six representative times (given above each panel in units of tg). Superposed on the 
color map of column density are velocity vectors (with normalization shown above each panel, where Cs is the isothermal sound speed), 
contours of critical flux-to- mass ratio, and locations of the stars formed (denoted by crosses). The number of stars A', is given above each 
panel. The simulation box has 10 Jeans lengths (Lj = 0.37 Tio/Ay pc) on each side. 
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Fig. 3. — Mass fraction of the magnetically supercritical material (upper curve) created in the initially subcriticaJ cloud of the standard 
simulation through ambipolar diffusion as a function of time. Only a fraction of this supercritical material resides in dense regions above 
ten times the average column density. The mass fraction of this dense, supercritical material is shown as the lower curve. 
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Fig. 4. — Snapshots of a variant of the standard simulation with a weaker outflow of / = O.Of (Model S2 in Table f). Superposed on the 
color map of column density are velocity vectors (with normalization shown above each panel, where Cs is the isothermal sound speed), 
contours of critical flux-to- mass ratio, and locations of the stars formed (denoted by crosses). The number of stars Nt is given above each 
panel. The simulation box has 10 Jeans lengths (Lj = 0.37 Tio/Ay pc) on each side. 
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Fig. 5. — Efficiencies of star formation of the six variants of the standard simulation hsted in Table 1. Tlie solid lines are for Models Rl 
and R2 that have the same power spectrum as the standard model, but different random realizations. The dotted lines are for Models II 
and 12 whose turbulent velocity fields are incompressible. The dashed (dash-dotted) line is for Model Nl (N2) that has a steeper (shallower) 
power spectrum than the standard model. 




Fig. 6. — Efficiencies of stcir formation for clouds of the same initial turbulent velocity field but different degrees of magnetization, 
measured by the dimensionless fiux-to-mass ratio To labeled beside the curves. The heavy (thin) solid lines are for strong (weak) outfiow 
cases with / = 0.1 (0.01). 
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Fig. 7. — Snapshots of a non-magnetic (Model Ol; panels a-c) and marginally supercritical cloud with Fq = 0.8 (Model Ul; panels d-f) at 
three representative times (given above each panel in units of tg). Superposed on the color map of column density are velocity vectors (with 
normalization shown above each panel, where Cs is the isothermal sound speed) and locations of the stars formed (denoted by crosses). 
The number of stars A'^, is given above each panel. The simulation box has 10 Jeans lengths (Lj = 0.37 Tio/Ay pc) on each side. 
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Fig. 8. — Efficiencies of star formation for the standard subcriticaJ cloud (Fq = 1.2; solid lines) and a magnetically critical cloud (Fq = 1; 
dotted lines) for different values of turbulent Mach number M. The role of strong turbulence in accelerating star formation is evident. 
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Fig. 9. — Spatial distribution of the dense cores identified in the standard simulation at the time t = 2 tg. The contours of column 
density are at S = 6, 10 and 14 Sq. Dense cores are enclosed by at least two contours. The core numbers are labeled. 
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Fig. 10. — Enlarged views of the dense cores shown in Fig. |5] except for panels (f) and (h), which display the core No. 6 at two different 
times (labeled above the panels). Plotted in each panel are the contours of column density and vectors of the velocity relative to the density 
peak (marked by a cross). The velocity normalization is shown above each panel, where Cb is the isothermal sound speed. 
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Fig. 11. — The mass- weighted rms Mach numbers of the gas (solid hne) and stars (dotted line) for the standard simulation as a function 
of time. Note the slowdown of turbulence decay at late times. 




Fig. 12. — The mass-weighted rms Mach numbers for the gas in the Fq = 1.2 (subcritical; solid lines), 0.8 (supercritical; dotted lines) 
and (non-magnetic; dashed lines) clouds as a function of time. The thick (thin) lines are for the strong (weak) outflow cases with 
/ = 0.1 (0.01). Evidently, appreciable replenishment of turbulence requires both a strong outflow and a high rate of star formation, a 
condition that is met only with the combination Fq = and / = 0.1. 



